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ABSTRACT 



Aims. The main sulfur-bearing molecules OCS, H 2 S, SO, S0 2 , and CS have been observed in four high mass dense cores (W43-MM1, 
IRAS 18264, IRAS 05358, and IRAS 18162). Our goal is to put some constraints on the relative evolutionary stage of these sources 
by comparing these observations with time-dependent chemical modeling. 

Methods. We used the chemical model Nahoon, which computes the gas-phase chemistry and gas-grain interactions of depletion and 
evaporation. Mixing of the different chemical compositions shells in a ID structure through protostellar envelope has been included 
since observed lines suggest nonthermal supersonic broadening. Observed radial profiles of the temperature and density are used to 
compute the chemistry as a function of time. 

Results. With our model, we undeiproduce CS by several orders of magnitude compared to the other S-bearing molecules, which 
seems to contradict observations, although some uncertainties in the CS abundance observed at high temperature remain. The 
OCS/SOt, SO/SO2, and H2S/SO2 abundance ratios could in theory be used to trace the age of these massive protostars since they 
show a strong dependence with time, but the sources are too close in age compared to the accuracy of chemical models and obser- 
vations. Our comparison between observations and modeling may, however, indicate that W43-MM1 could be chemically younger 
than the three other sources. Turbulent diffusivity through the protostellar envelopes has to be less efficient than 2 x 10 14 cm 2 s" 1 . 
Otherwise, it would have smoothed out the abundance profiles, and this would have been observed. 

Conclusions. The sulfur chemistry depends strongly on the ID physical conditions. Any observed set of abundances should be com- 
pared with a chemical model computed with the same temperature and density traced by the observations. In our case, no conclusion 
can be given on the relative age of IRAS 18264, IRAS 18162 and IRAS 05358 except that they are very close. W43-MM1 seems 
younger than the other sources. Turbulent mixing could occur in young high-mass protostars on a geometric scale that is too small to 
affect the chemical abundance profiles through the envelope. 

Key words. Astrochemistry - Turbulence - Stars: formation - ISM: individual objects: IRAS 18264, IRAS 18162, IRAS 05358 and 
W43-MM1 - ISM: molecules 



1. Introduction 

Despite our poor knowledge of the sulfur chemistry in the inter- 
stellar medium, the abundance of S-bearing molecules observed 



the sources being studied (Charnlev[ 1997 


Hatchell et all 1998 


Buckle & Fuller[ |2003|:| van der Tak et all 
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2004- |Herpin et all 120091) . Most of the rate coefficients of the re- 
actions involving sulfur-bearing molecules have not been stud- 
ied in the laboratory because of the difficulty of conducting ex- 
periments with them. Rate coefficients are then uncertain and 
many must be taken with caution. Another problem related to 
sulfur chemistry is the depletion of the element. Observations 
of the atoms in the diffuse medium have shown a relationship 
between the amount of depletion of the elements (in a more or 
less r efractory form) with the density of the cloud (see LfenkinsL 
2009, and references therein). Sulfur is one of the exceptions 
that do not show this depletion. As far as observations can go, 
the atomic abundance seems to be the same as the cosmic sulfur 
abundance of about 1 . 5 x 10~ 5 compared to the proton density 
(ISofia & Meveri. l200lh . If using such an elemental abundance 
in chemical models for cold clouds, we would overpredict gas 
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phase a bundances o f SO, SO2, and CS by several orders of mag- 
nitude dRuffle et all [1999). For this reason, a depletion of this 
element is usually assumed to proceed in dense clouds where 
the amount of S available for the chemistry is taken to be two 
orders of magnitude less than the cosmic abundance. Recently, 
by rean a lyzing hundreds of observations of atomic sulfur line, 
Jenkins (2009) has argued that not observing the sulfur deple- 
tion in diffuse medium is an observational bias. If confirmed, 
this result would solve one of the interstellar sulfur mysteries. 

The form of sulfur on interstellar grains is not known. 
OCS is the only mo lecule actually observed on interstellar ices 
dPalumbo et all[l997D . Since hydrogenationis expected to be the 
most efficient process on grains, H2S is also probably formed, 
although with a smaller abundance tha n the observational limit 
of 10 7 (Ivan Dishoeck & BlakeL Il998l) . H 2 SQ 4 has been pro - 
posed a nd looked for without a ny success (IScappini et all l2003). 
Finally, Wakel am et al.l |2004) propose that sulfur is in the form 
of polymers or aggregates of sulfur that would be quickly con- 
verted into atomic sulfur once evaporated in the gas phase. In 
protostellar envelopes, when the temperature increases, the sul- 
fured molecules evaporated fr om the gra i ns are t ransformed into 
SO fi rst and then into S0 2 (ICharnlevi Il997t IWakelam et all 
2004). Thus abundances of these species have been used on sev- 
eral occasions to trace the chemical evolution of such objects. 
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Ivan der Tak et al.l J2003) have studied the sulfur chemistry in 
nine high-mass protostars and find that high-energy transitions 
probing the inner parts of the protostars where the temperature 
exceeds 100 K were mandatory for using S -bearing molecules as 
chemical clocks. They were only able to obtain the abundance 
of SO2 in these regions. In addition, they propose that OCS is 
the main carrier of sulfur on grains, based on the high excita- 
tion temperature of the molecule and its high abundance in the 
protostars. 

In the present work, we made a ID modeling of the sul- 
fur chemistry in four high-mass compact objects, IRAS 18264, 
IRAS 18162, IRAS 05358, and W43-MM1, presumably rep- 
resentative of the youngest stages of massive star formation. 
We used the density and temperature profiles determined from 
a detailed modelin g of the SEP observ ed in these sources by 
iHerpin et al.l d2009h . lHerpin et all (120091) also observed the main 
S-bearing molecules in these sources in order to get an estima- 
tion of their abundance and possible evolutionary stage. The ob- 
served line profiles have shown a non-thermal line broadening 
that they attribute to supersonic diffusion. We then included a 
treatment of the diffusion in our chemical modeling in order 
to study its impact on the abundance profiles computed by the 
model. 

This paper is organized as follows. Our ID chemical model 
and the model parameters are described in section |2] Section [3] 
shows the abundance (compared to H) and abundance ratio pro- 
files computed by the model in the studied sources, as well as 
the effect of turbulent mixing. We discuss the particular prob- 
lem of CS in section |4] We compare our model predictions with 
the observed abundances in section and conclude in the last 
section. 



2. Model description 

The model u sed for this study is t he Nahoon chemical model 
developed by Wakelam et al. ( 20051 l20loh . This code originally 
computes the chemical evolution of a list of species as a function 
of time for a fixed temperature and density for a single point. In 
this new version of Nahoon, the chemical evolution is computed 
in a ID structure and accounts for diffusive mixing induced by 
turbulent transport. Starting from the chemical composition of a 
molecular cloud (see j32.31 l. we then compute the chemical abun- 
dances as a function of the radius to the center of the protostars 
and as a function of time using the physical structures (density 
and temperature profiles) observed in the four high-mass pro- 
tostars described in section 12.41 We compute the evolution of 
chemical abundances x, following 



dxj 



1 d 



d 



U + —f—Kpr —Xi 
pr^ or or 



(1) 



where f, and L, are the chemical production and loss terms, p is 
the density, k the turbulent diffusivity, and r the spherical radius. 
We compute turbulent diffusion in spherical coordinates. 



2.1. The chemistry 

Gas-phase chemistry is computed and species are allowed to 
deplete on grains and evaporate by direct temperature effect 
and indirect heating by collisions between grains and cos- 
mic ray particles. No reactions on grains are taken into ac- 
count in this work. Only an approximati on for the fo r matio n 
of H2 on grains is considered following iLeung et"ai1 (Q984). 



ey, 



The freezeout of ga s-phase species onto g rains is computed us- 
ing equation 1 of lHasegawa et al] (119921) . The evaporation of 
species from the grain mantles is also computed according to 
Hasegawa et al. ( 1992) for direct thermal evaporation and from 
lHasegawa & Herbstl (1 1993b for evaporation induced by colli- 
sions between grains and cosmic -ray particles. 

The chemi s try is described by the Srates network from 
Wakelam et al . (2004). Srates is a reduced network dedicated to 
sulfur chemistry in warm sources. It contains 929 chemical re- 
actions involving 76 gas-phase and 32 surface species for the 
elements H, He, C, O and S. This network has been assembled 
from different sources (OSU, UDFA and the NIST databases) 
and updated using the osu latest network (osu-09-2008, 
http://www.physics.ohio-state.edu/~eric/research.html). To se- 
lect the reduced net w ork, s everal previous work s wer e 
followed dRuffle et all 120021: iHollenbach & McKeel [197 
i Hartauist et all 1 19801: iPineau Pes Forets et all Il993b ICharnle 
1997) and the results of simulations obtained with this reduced 
network were tested against the ones of larger networks valid 
at low temperature. Srates can be downloaded at the following 
address: http://kida.obs.u-bordeauxl.fr/models. The binding en- 
ergies for OCS and H 2 S are taken to be 2888 and 2743 K, re- 
spectively. These bi nding energies were c omputed using equa- 
tion 1 described in Collings et all d2004l) and the evaporation 
temperatures of these species measured by TPP experiments 
from the same paper (Herma Cuppen, private communication). 
H2 cosmic -ray ionization rate is 1.3 x 10~ 17 s _1 . 



2.2. Diffusion 

Nonthermal broadening of lines is observed in various prestellar 
cores. This behavior is usually assumed to be the consequence 
of some sort of unresolved gas motions, whose exact origin and 
properties are unknown. However, we assume that these motions 
produce some kind of mixing inside the source and are often 
referred to as turbulence. Our purpose is to investigate the ex- 
act influence of a turbulent-like mixing on the distribution of S- 
bearing species, without presuming the exact nature of this "tur- 
bulence". The simplest w ay to do so is to use a mixing length 
approach (|PrandtlL fT925). We added mixing terms to all the dif- 
ferential equations in order to couple each spatial point in the 
protostar envelopes with its adjacent points. Turbulence induces 
an enhanced diffusion in the source, whose effective diffusivity 
can be written as 



(2) 



where vj and lj are a characteristic turbulent velocity and 
length scale (the mixing length), respectively. Expressing V7 and 
lj would require a detailed model for turbulence in protostel- 
lar cores. This is not, however, within the scope of the paper. 
Instead, we aim to use observed quantities as upper limits for 
both parameters. To estimate the maximum efficiency of the tur- 
bulent diffusion, we assume that the velocity Mj cannot be higher 

than the sound speed ( J~^, with R the gas constant, T the gas 

temperature, and [i the mean molecular weight), and the mixing 
length It is expected to be smaller than the radius of the enve- 
lope (~ 3 x 10 17 cm). We obtain a maximum turbulent diffusivity 
Kmax of about 2 x 10 22 cm 2 s _1 . 

2.3. Initial conditions for the protostellar envelope 

Protostars form from dense clouds. As a consequence, for the 
initial conditions of the protostars, we computed the chemi- 
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Table 1. Grain mantle composition. 



Species 


Abundance (/H) 


Ref. 


CO 


4 x 1CT 6 


(1) 


ii 2 u 




(2) 


H 2 CO 


2 x icr 6 


(3) 


CH 4 


5 x icr 7 


(4) 


CH3OH 


2 x 10~ 6 


(5) 


OCS 


5 x icr 8 


(6) 


H 2 S 


5 x 10~ 8 




S (mod 1) 


1.45 x irr 5 




S (mod 2) 








References: (1) 



Gibb et alj <200Cb. (2 ) ISchutte et alj 



iKeane et alJteOOll). (4)IChiar et al]t 19961) 
(6) Palum bo etal.U1997h 



Booeert et al. 



1996), (3) 
(1998:) and 



cal composition of a dense molecular cloud with a temperature 
of 10 K, a total H density of 2 x 10 4 cm 3 , a visual extinc- 
tion of 30, and a cosmic -ray ionization rate of 1.3 x 10~ 17 s _1 . 
For the dense cloud, species are assumed to be initially in the 
atomic form (coming from the diffuse medium) except for H2. 
Elemental abundances (compared to total H) used for this calcu- 
lation are 0.09 for He, 2.56 X 10" 4 for O and 1.2 x 10~ 4 for 
C + dWakelam & Herbsll 12008b . For S + , we chose an elemen- 
tal abundance of 5 x 10~ 7 for the dense cloud in order to ap- 
proximately reproduce th e observed abundance s of gas-phase 
S-bearing molecules (see IWakelam et all 12004]) . Using our re- 
duced network, we have to integrate over a long time (10 7 yr) to 
obtain abundances that are similar to the ones observed in dense 
clouds . We in fact obtain ab undances similar to our previous 
study (IWakelam et all 12004 Composition A of Table 1). This 
time is probably too long compared to the age of dense molec- 
ular clouds. Using larger networks of more than four thousand 
gas-phase reactions, observations in dens e clouds are usually re - 
produced at a few 10 5 yr (see for instance Wakela m et al.l l2006). 
The reason we need more time is most likely that we have fewer 
chemical species and that it takes more time for the ionization 
fraction of the gas to decrease. At high temperatures however, 
we do not have this problem. This calculation gives us our gas- 
phase composition prior to the formation of the protostar. The 
CO abundance computed by the model on the grain surfaces (re- 
sult of the depletion of gas-phase CO) is about ten times larger 
than what is observed. We decreased this value to 4 x 10" 6 (CO 
abundance observed in inte rstellar ices towa rds the massive pro- 
tostar line of sight W33 bv lGibb et aLT 120001) and put the rest of 
the CO back in the gas phase. The total CO gas-phase abundance 
is then 6.08 x 10~ 5 . This change will not affect our chemistry 
since most of the protostar envelopes are above the evaporation 
temperature of CO (~18 K). Like most dense cloud chemical 
models, our simulations fail to reproduce the observed limits on 
the O2 abundanc e in dense clouds, which is 5 x 10 s (compared 
to H) in L134N dPagani et alll2003l) . 160 times smaller than the 
abundance computed by our model. We then decreased the gas- 
phase abundance to its observed limit, and the rest of the oxygen 
we put in the abundance of atomic oxygen. The computed O2 
abundance on grains is 5 x 10~ 5 . For H2O, our model predicts a 
gas-phase abundance of 4 x 10" 8 , below the observational limit 
in L134N (1.5 x 10' 7 ISnell et"aTl l2000t) so we did not change it. 

Since we did not compute the grain surface chemistry occur- 
ring during the dense cloud phase, we modified the grain surface 
composition obtained from our modeling. For the species ob- 
served on inters tellar grains, we adopt ed the abundances listed 
in Table [TJ (see Wakela m et all [2004). Note that C0 2 , one of 
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Fig. 1 . Temperature and density profiles in the four sources from 
HerpinetalJ(|2009). 



Table 2. S-bearing molecular abundan ces (xlO 
the four sources bv lHerpin et al.l (120091) . 



■10 



) derived in 



Molecule 


T(K) 


w43-mml 


Sources 
IRAS 18264 IRAS 05358 


IRAS 181f 


CS 


60 


36 


4.0 


2.3 






100° 


160 


1.1 




4.6 


OCS 


60 


92 


2.5 


2.4 


1.4 




100 


130 


1.1 


3.7 


3.3 


H 2 S 


60 


3.0 


1.3 


1.2 


3.1 




100 


4.0 


0.15 


0.51 


0.30 


SO 


60 


0.90 


0.48 


1.1 


2.4 




100 


24 


0.64 




3.6 


S0 2 


60 


2.0 


0.42 


0.62 


1.7 




100 


160 


2.0 


6.4 


7.6 



" The CS abundance at this excitation temperature was derived usin; 
the C S (7-6) intensity line blended with an H 2 CO line (see Herpin et al 
I2009h . 

The symbol - means that the no abundance was derived because the 
molecular line has not been observed. 



the main constituent of grain mantles, is not present in our 
network. OCS is the only S-bearing species observed in the 
solid state in the interstellar medium. H2S is highly suspected 
of being present, although at a nondetectable level. We set its 
abundance on grain surfaces to the upper limit calculated by 
Ivan Dishoeck & Blakd (119981) . All the rest of the sulfur was as- 
sumed to be either in the atomic form on grains, our model 1, or 
in a refractory form that cannot evaporate at the temperatures of 
hot cores, our model 2. This gas-phase and grain mantle compo- 
sition is used as initial conditions for the whole envelope of the 
four protostars we studied. Since we do not include any chem- 
ical reaction on the grain surfaces, the gradual depletion of the 
species during the evolution of the chemistry only affects the 
fraction of species in the gas and on the surface and does not 
change the chemistry itself. 

2.4. Selected sources: physical structure and observed 
abundances 

The source sample from Herpin et al. (2009) is made of two 
IR-quiet dense core s (MSXFlux7i„ m < 10 Jy at 1.7 kpc, follow- 
ing the definition of Mot te et all 120071) and two slightly brighter 
ones (IRAS05358+3543 and IRAS 18 162-2048), with bolomet- 
ric luminosities (0.6-2.3)xl0 4 L at distances 1.8-5.5 kpc and 
sizes of ~0.1 1—0. 13 pc. These massive dense cores are meant to 
represent the earliest phases of the high-mass star formation. 

From the fitted SEDs and the literature, Herpin et al (2009) 
propose a rough evolutionary classification of the four ob- 
jects (W43-MM1 being the youngest object): W43-MM1 -> 
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IRAS 18264- 1152 -» IRAS05358+3543 -> IRAS18162-2048. 
Nevertheless and despite this youth, W43-MM1 appears to lie 
apart from the other massive dense cores as it has very likely 
already developed a hot core. 

The temperature and density profiles for the four sources 
have been derived by iHerpin et al.l (120091) from the analysis of 
the spectral energy distribution (S EP) with the radiative trans- 
fer co de MC3D (IWolf et all 1 19991) (see also iMarseiUe etafl 
2008). The profiles are shown in Fig. [TJ Temperatures and 
densities increase towards the center of the sources. The 
three IRAS sources show very similar profiles, whereas W43- 
MM1 is denser and colder. For the chemical modeling, we 
assumed that the gas temperature and density are abruptly 
raised from the cold molecular phas e to the ones currently 
observed following previous wor ks dDotv & Neufeldl [1997; 
Ivan der Tak & van DishoeckT, 120001; iDotv et all 120021) . This is of 
course a strong assumption but we prefer not to add another un- 
certain parameter. Considering the fast evolution of sulfur chem- 
istry at high temperature (see section 3.1), we do not expect this 
assumption to affect our main conclusions. Only a careful cou- 
pling of the chemistry with a dynamical infalling envelope can, 
ho wever, answer that question. 

IHerpin et all (120091) observed several transitions of OCS, 
H2S, SO, SO2, and CS in the four sources, as well as some 
of their 34 S minor isotopologues. They constrained the species 
abundances at different radii depending on their excitation tem- 
peratures (60, 75, and 100 K) (see Table 5 of their paper). We 
report in Table [2] their observed abundances for excitation tem- 
peratures of 60 and 100 K. We assume that we are at LTE so 
the excitation temperature of the gas reflects the gas and dust 
temperatures. The presence of a hot core in W43-MM1 might 
explain the observed higher abundance of H2S and OCS in that 
source relative to the other ones. Moreover, strong shocks within 
this source probably influence the chemistry. One surprising re- 
sult is the small abundance of the S-bearing molecules found 
in all four sources compared to the cosmic abundance (three to 
four orders of magnitude smaller). We compare our modeling 
with their observational results below. 



3. Modeling results 

3.1. Abundances through the protostellar envelopes 

Figures [2] and [3] show the abundances of the main S-bearing 
molecules computed by the chemical model as a function of 
the radius to the central object in the envelope of the four pro- 
tostars at 10 4 yr which was approximately the age derived by 
IDotv et alJ d2002l) for the massive protostar AFGL 2591. This 
age was derived by comparing the predictions of chemical mod- 
els with a large number of molecules observed in this source. 
Although AFGL2591 may be a bit more evolved than the four 
sources studied here (AFGL2591 is a mid-IR bright HMPO, 
while the 4 others are mid-IR quiet HMPOs), all these are still 
HMPOs, i.e. not hot molecular cores or ultra-compact HII re- 
gions. As far as we know, it is the only large chemical study of 
such sources so we will use it as a reference. The two figures 
have been obtained for a different amount of the atomic sulfur 
initially available on the grain surfaces (models 1 and 2, see !2.3l ). 
IRAS 18264, IRAS 18162 and IRAS 05358 show very similar 
molecular profiles since their temperature and density gradients 
are similar (Fig.[TJ, except for radii larger than 7 x 10 3 AU where 
the depletion of the species is more pronounced in IRAS 18162. 
W43-mml shows different abundance profiles, especially the 
outer envelope where the molecules are more depleted because 
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Fig. 2. OCS, H 2 S, CS, SO and S0 2 abundance (/H) profiles as a 
function of radius from the center in the envelope of four high 
mass protostars (IRAS 18264, IRAS 18162, IRAS 05358 and 
W43-MM1). Initial composition of Model 1 is used. Integrated 
time is 10 4 yr. 
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Fig. 4. OCS, H 2 S, CS, SO and S0 2 abundance (/H) profiles as 
a function of radius from the center in the envelope of IRAS 
18264. Initial composition of Model 2 is used. Integrated times 
are 10 3 yr for the figure on the top and 5 x 10 4 yr for the figure 
on the bottom. 
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Fig. 5. Abundance (JH) profiles of a selection of species as a 
function of radius from the center in the envelope of W43-MM1 . 
Initial composition of Model 2 is used. Integrated time is 10 4 yr. 



of higher densities. In this source, both models 1 and 2 lead to 
smaller abundances of SO, H 2 S, and CS molecules compared 
to the other sources. All species show strong variations in their 
abundances with radius. 

The radius at which the gas and dust temperatures are around 
40-50 K represents a transition from an "outer envelope", where 
the S-bearing molecules are depleted on the grain mantles, and 
an "inner envelope" where the abundances of the S-bearing 
species are much larger. In the inner envelope, towards the cen- 



Fig.3. OCS, H 2 S, CS, SO and S0 2 abundance (JH) profiles as a 
function of radius from the center in the envelope of four high 
mass protostars (IRAS 18264, IRAS 18162, IRAS 05358 and 
W43-MM1). Initial composition of Model 2 is used. Integrated 
time is 10 4 yr. 
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Fig. 6. Abundance (/H) profiles of a selection of species as a 
function of radius from the center in the envelope of IRAS 
18264. Initial composition of Model 2 is used. Integrated time 
is 10 4 yr. 

ter of the protostar, all molecules except SO2 again show a de- 
crease in their abundances at radii depending on the source, the 
molecule and the model. In IRAS 18264 and model 1 for in- 
stance, OCS drops at 200 AU and SO abundance at 400 AU. 
The H2S abundance decreases by more than one order of mag- 
nitude between 1000 and 400 AU before increasing again. In all 
sources, SO2 is the only species whose abundance is rather con- 
stant in the inner envelope, whereas the CS abundance is quite 
low and shows a strong dip where the temperature is about 200 K 
(see discussion in section @}. In the very outer part of the en- 
velopes (at radii larger than 10 000 AU), the species abundances 
increase again because the density of the gas is lower and species 
have not had time yet to deplete on the grains. 

In model 1, a large quantity of atomic sulfur is evaporated 
in the gas phase when the temperature is above ~20 K. Atomic 
sulfur is then quickly converted into SO by reacting with O2. SO 
in turn produces SO2 by reactions with OH and O. For temper- 
atures higher than ~ 40 - 50 K, OCS and H2S are evaporated 
in the gas phase, and their abundances are determined by their 
initial solid abundances. In model 2, it is the destruction of H2S 
through the following paths that releases the necessary sulfur 
to produce SO and SO2 at high temperature: H2S + H30 + — > 
H 3 S+ + H 2 0, H 3 S + + e " -» HS + H 2 and HS + O — > SO + H. 
SO2 is the dominant S-bearing molecule after a few 10 4 yr in all 
models at radii smaller than 7 x 10 3 AU. We anticipate the com- 
parison with the observation given in section [5] by saying that 
model 1 is far from reproducing our observations since very low 
abundances of S-bearing molecules have been observed. For this 
reason, we only discuss model 2 in the rest of the paper. 

Some other species abundances (for model 2 and 10 4 yr) are 
shown in Fig.[5]for the source w43-mml and Fig.|6]for the source 
IRAS 18264 (representative of the IRAS sources of our sample). 
In IRAS 18264, the CO abundance is constant over the envelope 
(~ 8 x 10~ 5 ) since the temperature is above 20 K, CO deple- 
tion occurs whereas at radii larger than 10000 AU in w43-mml. 
Inside the protostellar envelopes (radii smaller than 200 AU for 
w43-mml and 600 AU for IRAS 18264), the main reservoir of 
oxygen is water, except for CO. At larger radii, the oxygen is 
mainly in the atomic form. This transition between H2O and O 
corresponds to the dip in the CS abundance (see also section 4). 
Similarly, the peak in HOCS + , precursor of CS, around 2000 AU 
in w43-mml and 5000 AU in IRAS 18264 is responsible for the 
peak in the CS abundance. 

The abundances in Figs. |2]and[3]are displayed for an age of 
10 4 yr. By age, we mean the time during which the chemistry 
evolves after the temperature and density of the envelope have 
reached their present state. As an example, the chemical pro- 
files in the source IRAS 18264, for the elemental abundances of 
model 2, are shown in Fig. [4] for two different chemical ages: 




Fig. 7. SO/SO2, OCS/SO2, H2S/SO2 and CS/SO2 abundance ra- 
tios as a function of time at radius where the temperature is 
100 K in IRAS 18264 and W43-MM1. Initial composition of 
Model 2 is used. 



10 3 and 5 x 10 4 yr. The chemistry evolves quickly inside the 
protostellar envelopes and the abundances are similar at 10 3 and 

10 4 yr except for the SO2 abundance, which increases signifi- 
cantly with time. As time evolves after 10 4 yr, OCS and CS at 
radii 2000-6000 AU are destroyed by HCO + . Abundances in the 
outer parts of the envelope also evolve toward more depletion. 

It is interesting to notice that the abundance of OCS does not 
change much at radii smaller than 2000 AU with time. In fact, 
it takes more than 10 s yr to significantly destroy this molecule 
in the conditions of the studied massive protostars. If OCS were 
th e major carrier of S-beari ng molecules on grains as suggested 
by Ivan der Tak et all (120031) . time scales for the sulfur chemistry 
would be longer than what we considered here. 

3.2. Abundances ratios 

It is common when comparing with observations to use abun- 
dance ratios between two observed species rather than abun- 
dances compared to H2. We have observational constraints on 
the abundances at two different radii corresponding to gas tem- 
peratures of 60 K and 100 K (see section [2~4b . The abundance 
(and abundance ratios) observed at 60 K are probably not good 
tracers for the evolution of the protostars since we probe the re- 
gions of evaporation versus depletion processes of the sulfur- 
bearing species. To simplify the problem since we have many 
parameters, we display in Fig. [7] the abundance ratios of SO, 
OCS, H 2 S, and CS compared to S0 2 in IRAS 18264, and W43- 
MM1 as a function of time for the radii corresponding to the 
temperature 100 K (2680 AU in IRAS 18264 and 640 AU in 
W43-MM1). We did not show the abundance ratios for IRAS 
18162 and IRAS 05358, because they are very similar to those 
in IRAS 18264. All four abundance ratios decrease with time. 
SO/SO2 shows a monotone variation with time contrary to the 
three other ratios, which drastically decrease at times between 
10 4 and 10 5 yr. Since we do not reproduce the observed abun- 
dance of CS, OCS/SO2 and H 2 S/S0 2 seem to be the best candi- 
dates to constrain some evolutionary time for these sources. 

3.3. Influence of diffusion 
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Fig. 8. OCS, H 2 S, CS, SO and S0 2 abundance (/H) profiles as 
a function of radius from the center in the envelope of IRAS 
18264. Initial composition of Model 2 is used. Integrated times 
is 10 4 yr. Three different turbulent diffusivity k are used for each 
panel: 2 x 10 14 crrT 2 s _I for the upper panel, 2 x 10 16 cm~ 2 s _I 
for the panel in the middle and 2 x 10 18 cirT 2 s _1 for the lower 
panel. 



The general effect of diffusion is to smooth out the abun- 
dance profiles across the envelope and decrease the depletion of 
the species abundances in the outer parts. Abundance profiles for 
the source IRAS 18264 and model 2 are shown in Fig.|8]for three 
different values of the turbulent diffusivity k: 2 x 10 14 cirT 2 s _1 , 
2 x 10 16 cirT 2 s _1 , and 2 x 10 18 cirT 2 s _1 . These figures can be 



compared to Fig. [3] where k is zero. Our largest k is four orders 
of magnitude smaller than the maximum value derived in sec- 
tion l2.2l but the effect on the chemistry is already strong. In that 
case, the abundances of OCS, H 2 S, SO, and S0 2 are then con- 
stant and equal up to radius 5 x 10 3 AU. The CS abundance is 
also rather smooth, but four orders of magnitude less than the 
other species. For smaller k, the smoothing is less efficient and 
becomes negligible for k lower than 2 x 10 14 cm~ 2 s _1 . 

4. The CS abundance 

As already noticed, we predict that CS is much less abundant 
than the other S-bearing species contrary to the observations, 
whic h show abundances o f the same order for all species (see 
also Wakelam et a In our model, CS is produced at 

high temperature by the dissociative recombination of HCS + 
in the central part where the temperature is above 300 K and 
the dissociative recombination of HOCS + at lower temperature. 
HCS + is itself produced by C + + H 2 S and HOCS+ by HCO + 
+ OCS and + OCS. For the destruction, we also have two 
regimes: CS + H + at temperatures hi gher th an 300 K and CS 
+ O fo r lower temperatures. ICharnlevI (ll 9971) and lHatchell et al .1 
( 1998) have produced much more CS than in this work because 
they did not include any atomic oxygen in their initial condi- 
tions. Except for helium and hydrogen, they considered that all 
species were already in a molecular form. Oxygen was then in 
the form of CO, H 2 Q , and 2 . Contrary to what was said in 
IWakelam et all d2004t) . changing the rate coefficient of the dis- 
sociation of CO by secondary photons produced by cosmic -rays 
does not change the abundance of CS predicted at high temper- 
ature. Sinc e this problem is not s een in cold dense sources (see 
for instance lWakelam et al.Ll2006h . we may be missing some for- 
mation paths for CS at high temperature. Some rate coefficients 
may also not be accurate for the temperature range where we 
are using them. As an example, some new recommendations for 
the rate coefficient of the neutral-neutral reaction O + CS for 
temperatures between 150 and 300 K have been posted on the 
KIDA database (http://kida.obs.u-bordeauxl.fri following this 
work. However, the small proposed changes in this rate coef- 
ficient does not change the gas-phase CS abundance in our sim- 
ulations. 

That we do not reproduce the CS abundance at 100 K may 
also come from our other assumptions in our modeling such as 
our initial conditions, the temperature and density profiles, the 
brutal rise in the temperature and density from the molecular 
cloud stage. This, however, should not influence the abundance 
of the other S-bearing molecules by more than a factor of a few 
since they all have similar observed abundances. The CS prob- 
lem may also be related to the form of sulfur on grains. Recent 
laboratory experiments show that H 2 S on ISM analog surfaces 
are easily destroyed by energetic particles, and surface S then re- 
combines to form OCS, SQ 2 , C 2 S, and a majority of sulfur-rich 
residuum dGarozzo et al.L l2010h . The surface C 2 S may then be 
evaporated in the gas phase and contribute to the CS chemistry. 

Last but not least, we stress that the discrepancy might also 
originate in the observed abundances. The observed abundance 
of CS representative of the 100 K layers in the protostellar en- 
velopes was derived using the CS (7-6) line emission, which 
has an upper level energy of 65.8 K. First of all, these abun- 
dances are derived independently for each observed transition of 
each species, assuming that the emission for a given line mostly 
comes from one single place at a given temperature and den- 
sity, whereas lines are excited over several layers. Considering 
the upper level energy of CS (7-6), it is likely that this transi- 
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tion is also excited in the layers of the protostellar envelopes at 
temperatures below 100 K where our chemical model predicts 
larger abundances of CS. The CS abundance is predicted by our 
chemical models to vary much more with the radius than for the 
other S-bearing abundances at these temperatures. As a result, 
an uncertainty about the location of the emission may lead to a 
larger uncertainty concerning the abundance. Furthermore, this 
line was blended with an H2CO line, and the procedure to sepa- 
rate lines did not w ork well in IRAS 18264 and w43-mml (see 
IHerpin et all 12009). The CS abundances at 100 K derived by 
( Herpin etall T2009) may then be overestimated, and one would 
need to observe higher energy transitions of CS (unfortunately 
not accessible from the ground) to confirm its abundance at high 
temperature. 

5. Comparison with observations 



In this section, we discuss the comparison between our mod- 
eling and observation s of S-bearing mole cules in the four high- 
mass protostars from iHerpin et al.l (12009b . As a summary, very 
small abundances have been found compared to the cosmic 
abundance of sulfur (1.5 x 10~ 5 compared to H). The sum of 
all observed S-bearing species is only 2 x 10~ 9 to 5 x 10~ 8 . This 
result suggests that we are closer to model 2 than model 1, i.e. a 
large part of the sulfur in this object is in a stable, not observed 
form and does not participate in the gas phase chemistry. We dis- 
cuss this point in the next section. We are far from reproducing 
all the observed abundances correctly. We overproduce them by 
more than one order of magnitude even using model 2, except 
for CS which is strongly underestimated. 

The OCS/SO2, H2S/SO2, and SO/SO2 abundance ratios 
could in theory be used to trace the age of these sources. We 
report in Fig. [9] the abundance ratios predicted by model 2 as 
a function of time for the radius corresponding to 100 K. As 
stated in section 3.2, the abundances derived at radii where the 
temperature is about 60 K trace regions of the envelopes where 
the depletion processes compete with the evaporation from the 
grains. For this reason we focus on the abundances obtained for 
shells of gas at higher temperature probably more representa- 
tive of the evolution of the protostar itself. In the same figures 
we plot the ratios observed in each source. The OCS/SO2 and 
H2S/SO2 abundance ratios are above the observed ones until a 
few 10 4 yr and then decrease sharply to the observational values 
and below. As a consequence, in the observed range for these ra- 
tios, the constraint on the age is very strong. This time depends 
on the physical properties of the source. The three IRAS sources 
being very similar, the time of decrease is almost the same. W43- 
MM1 is denser and the time at which the abundance ratios drop 
is less. The constraint on the relative age of the sources pro- 
vided by the observed ratios OCS/SO2 and H2S/SO2 is in good 
agreement with the evolution sequence W43-MM1 — » IRAS 
18264 -> IR AS 05358 -> IRAS 18162 proposed using other 
methods (see IHerpin etaTl l2009t). If one consid ers the uncer- 
tainties in both the observed (IHerpin et all [2009) and modeled 
( Wak elam et all 120051) abundance ratios, no definitive conclu- 
sion can be drawn from these observations. The constraints on 
W43-MM1 are maybe more robust and indicate that the source 
is younger than the three other ones. One important conclusion 
here is that more than the age of the source, it is the physical 
conditions in the envelope that directly determine the variation 
in the OCS, H2S, and SO2 abundance. There is of course a rela- 
tion between the age and the physical conditions. OCS/SO2 and 



H2S/SO2 can, however, not be used to compare protostars that 
are close in age. The SO/SO2 abundance ratios are more sensi- 
tive to time. We do need a very young chemical model for W43- 
MM1 to reproduce the observed abundance ratio in this source, 
which does not agree with the other constraints. The SO/SO2 
abundance ratios observed in IRAS 18162 and IRAS 18264 are 
also too close to conclude anything about their age. 

Observations also seem to indicate significant variations in 
the species abundances with the radius in a source. H2S, for in- 
stance, was observed to vary by one order of magnitude in IRAS 
18264 for radii between 2.6 x 10 3 and 5.4 x 10 3 AU, which 
is approximately reproduced by our model. Including diffusion 
smooths the predicted abundance profiles. If chemical variations 
are so strong across these objects, then we would expect the tur- 
bulent diffusivity to be smaller than 2 x 10 14 cm 2 s _1 . Supersonic 
velocities due to diffusion have been observed in these sources. 
Assuming a turbulent velocity of about 1 km s as observed, the 
maximum mixing length would be a few 10~ 4 AU, much smaller 
than what we would expect. The very small mixing efficiency 
associated with the observed nonthermal velocity raises ques- 
tions about the nature of this "turbulence". Either turbulence has 
a very peculiar nature, or more likely, the observed structure is 
unresolved and composed of smaller structures of varying turbu- 
lence intensity. Of course, nonthermal motions can be produced 
by a variety of processes (rotation, infall, outflow, etc), which 
will not always induce an efficient mixing, but can still be mis- 
taken for turbulence, at least partly. 

6. Conclusions 

We performed a ID chemical modeling of the sulfur-bearing 
species in four young high-mass protostars (IRAS 18264, IRAS 
18162, IRAS 05358, and W43-MM1) using tem perature and 
density profiles determined from their SED by IHerpin et al.l 
d2009l) . In addition to the time-dependent chemistry, we studied 
the effect of a turbulent-like mixing on the abundance profiles 
through the envelopes. Here are our main conclusions: 

- Sulfur chemistry depends strongly on the ID physical condi- 
tions. Any observed set of abundances should be compared 
with a chemical model computed with the same temperature 
and density traced by the observations. 

- To use sulfur chemistry as chemical clocks, observations 
tracing the gas at temperatures higher than 70 K should 
be used to avoid confusion with depletion mechanisms. At 
lower temperatures, S-bearing molecules stick on grains so 
that their abundance in the gas phase will also depend upon 
this process. 

- In our case, no conclusion can be drawn on the relative age 
of IRAS 18264, IRAS 18162 and IRAS 05358 if one con- 
siders the uncertainties in the observed and modeled abun- 
dances because these sources are too close in age. W43- 
MM1 seems, in contrast, younger than the other sources. 

- Turbulent mixing could occur in young high-mass protostars 
on a too small geometric scale to affect the chemical abun- 
dance profiles through the envelope. Such a small scale sug- 
gests that either turbulence is very unusual or the structure 
is unresolved. There is a possibility that at least some of the 
nonthermal broadening of lines would be the consequence 
of a non turbulent process, inducing a very inefficient mix- 
ing inside the protostellar envelopes. 

- The CS molecule is predicted by our chemical models to 
be less abundant than the other S-bearing species by sev- 
eral orders of magnitude at temperatures above 100 K (see 
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also IWakelam et all 120041) . This seems to disagree with 
the observat i ons, a lthough the observed abundances from 
iHerpin et all (2009) may be overestimated. If the disagree- 
ment were true, however, this suggests that some formation 
paths are missing in our networks at high temperature. 
- Our comparisons between observations and chemical mod- 
eling seem to indicate that the majority of the sulfur is 
still in a refractory form on grains. It is interesting to note 
that [Wakelam et al. (2004) also need a depletion by a fac- 
tor of ten of the elemental sulfur abundance to reproduce 
observations in Orion-KL, whereas no depletion is manda- 
tory to reproduce the observations in the low-mass proto- 
star IRAS 16293-24 22. Based on laboratory experiments, 
iGarozzo et all (1201 Oh have recently suggested that H2S on 
interstellar grains are easily destroyed by energetic particles 
to form OCS, SO2, and a majority of sulfur-rich residuum, 
which could be polymers of sulfur or amorphous aggregates 
of sulfur, as suggested bv lWakelam et al.l (12004). In the en- 
vironments of high-mass star formation, interstellar ices are 
probably exposed to stronger particle fluxes so that a larger 
quantity of the atomic sulfur could be converted in refractory 
form. The identificat ion of the sulfur-rich residuum found by 
IGarozzo et all d2010 ) could certainly bring new insight into 
the reservoir of sulfur in the interstellar medium. 
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Fig. 9. OCS/S0 2 (upper panel), H 2 S/S0 2 (middle panel) and 
SO/SO2 (lower panel) abundance ratios as a function of time 
at radius where the temperature is 100 K in the four studied 
sources. Initial composition of Model 2 is used. For each source, 
the observed abundance ratio is reported on the curves as hori- 
zontal thick line. Two values are drawn for IRAS 18162 as two 
SO2 abundances have b een determined using two different lines 
(see|Herpin et al., 2009). No observed value SO/SO2 is reported 
for IRAS 05358. 



